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We present the recent Fermilab calculations of the masses of the light quarks, using tadpole-improved 
Sheikholeslami-Wohlert (SW) quarks. Various sources of systematic errors are studied. Our final result for 
the average light quark mass in the quenched approximation evaluated in the MS scheme is rnq{iJ, = 2 GeV; Uf = 
0) = (m„ -f- md)/2 = 3.6 ± 0.6 MeV. 



1. Introduction 

We present recent results on the light quark 
mass determination using the SW action [1], 
which are updates of the last year's results [2]. 
For results from wilson and staggered fermions 
see [3] [4] [5] . 

The basic procedure is to extract the pseu- 
doscalar masses (mps) numerically for a range of 
quark masses and determine the linear coefficient 
in the chiral extrapolation. 
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Lattice details (n/ = 0) 
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where rrv^ = ln(l + l/2k— 1/2kc), with k = kuq 
and uq = \/<Up >mc [6]. 

Using the experimentally measured pion mass 
as an input, we obtain the light quark bare mass 
m[^"*, which is the average of the up and down 
quark masses. We convert it to the light quark 
mass TOqin the MS scheme by perturbation the- 
ory. 

Table 1 shows the lattices used for the sim- 
ulation. We use the SW fermion action. For 
/3 = 5.5, 5.7 and 5.9 the clover coefficient c is the 

'Presented by T.Onogi. 



tadpole improved tree-level value I/mq- However, 
for /? = 6.1, we use c= 1.40 instead of 1/mq = 1.46. 
All calculations are done in the quenched approx- 
imation. The lattice spacing a is determined from 
the IP-IS charmonium splitting. 

2. Systematic Errors 

We use the multi-state smearing method [7] 
to suppress excited state contamination. The 
smearing sources are fits to the measured wave- 
functions of the pseudoscalar ground and excited 
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Figure 1. The pseudoscalar effective mass plot 

for /3 = 5.9, K = 0.1385. The 1 x 1 one-state (dia- 
mond), 2x2 two-state (square), 3x3 three-state 
(circle) fits are shown. 



Figure 2. The Chiral extrapolation of m%g for 
/? = 5.9. The four k fit (solid line) and the three 
K fit (dashed line) are almost indistinguishable. 



states with the following forms. 



/is(r) = 
/2s(r) = 



exp(-/iisr), 

(1 - z/2s:)exp(-/x2sr). 



(2) 
(3) 



For (3 = 5.5, 5.7 and 5.9, wc use IS, 2S and lo- 
cal sources, while for (3 = 6.1, only IS and local 
sources are used. We choose 2x2 two-state fits 
as our best fits. In order to estimate the sys- 
tematic error of excited state contamination, we 
compare our best fits with the results from 1x1 
one-state and 3x3 three-state fits. We find that 
the difference is less than 1% for (3 = 5.7, 5.9 and 
about 1-1.5% for (3 = 6.1. (See Figure 1.) 

As the chiral extrapolation error, wc take the 
difference in the chiral extrapolation with three 
k's and four k's. The results are again less than 
1% for (3 = 5.7 and 5.9, and about 3% for (3 = 6.1. 
(See Figure 2.) 

The one loop the renormalization factor which 
connects the lattice bare mass with MS mass is, 



m,(/x) 



lat 



l+av(g*)(7o(lnC„-ln(a/i))) .(4) 



The mean-field improved bare mass m^"* is given 



by m[,"* = mo/(l 
theory [6]. 70 = 



-(7r/3)av+- • • ) in perturbation 
2/7r is the leading quark mass 



Table 2 

Our results for A, the coefficient of the lin- 
ear fit in the chiral extrapolation, mg"* (in 
MeV), the tadpole improved lattice bare mass, 
and mg(2GeV) (in MeV), in the WS scheme, 
renormalized at 2 GeV {q* = w/a, 1/a.) . 





5.5 


5.7 


5.9 


6.1 


A 

f^lat 

TOq(7r/a) 

TOg(l/a) 


5.3(7) 
4.34(17) 
4.75(19) 
6.20(25) 


4.1(5) 
3.9(1) 
4.41(12) 
4.89(13) 


3.05(7) 
3.3(1) 
3.90(13) 
4.19(14) 


2.24(10) 

3.2(1) 
3.84(18) 
4.05(19) 



anomalous dimension. Cm for SW-improved light 
quarks is 4.72 [8]. 

Using Eq.(4), we first convert the lattice quark 
mass to the MS mass at /z = tt/o or 1/a, then 
run it to the common scale of 2 GcV. In Eq.(4), 
there is another scale q* , which is the scale for the 
gauge coupling constant. Since we do not know 
the two-loop correction, it is not obvious which 
scale we should take for q*. We estimate the size 
of unknown higher order corrections to Eq. (4) by 
varying q* between 1/a and n/a. This procedure 
is consistent with assuming a coefficient of order 
unity for the ay term. Our estimates are 30%, 
13%, 7%, 5% for (3 = 5.5, 5.7, 5.9, 6.1. 
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Figure 3. The continuum extrapolation of Wiq. 

The upper value is the /? = 6.1 result, and the 
lower value is the naive linear extrapolation of the 
13 = 5.7, 5.9, 6.1 data. The data for /j, = q*=l/a 
(circle, solid line), 7r/a (diamond, dashed line) are 
presented in the same graph. 



There are both 0{a.a,) and O(a^) corrections 
to the action, and the continuum extrapolation 
could change depending on the relative size of 
these subleading terms. All we can say is that 
there is a systematic downward trend as we ap- 
proach to the continuum. Without a theoretical 
argument to tell us about the a-dependence, we 
take the /? = 6.1 result as an upper value and take 
the linearly extrapolated value using (3 = 5.7, 5.9, 
6.1 as a lower value. Our estimate of the contin- 
uum extrapolation error is 11%. (See Figure 3.) 



3. Summary 

In summary, our error estimates are, 

excited states < 1.5% 

chiral extrapolation ~ 3% 

perturbative 5% 

continuum extrapolation 11% 



combined 



17% 



The perturbative and a dependent errors are in- 
tertwined. We combine them linearly in the fol- 
lowing way. As we saw earlier, the scale of the 



coupling constant q* is arbitrary. When we dis- 
cuss the continuum limit, we therefore perform 
the extrapolation of the data for both q* = 1/a 
and TV /a (Figure 3). The outer points so obtained 
are taken as the limits of the combined error bar. 
The remaining errors are much smaller and com- 
bined in quadrature. Our final result for the light 
quark mass in the MS scheme in the quenched 
approximation is 



mgin = 2GeV;n/ = 0) = 3.6±0.6MeV. 
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